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Abstract 

We discuss the separability of the position-dependent mass Hamilto- 
nian in cylindrical coordinates in the framework of a radial power-law 
position-dependent mass. We consider two particular radial mass set- 
tings; a harmonic oscillator type, and a Coulombic type. We subject 
the radial harmonic oscillator type mass to a radial harmonic oscillator 
potential and the radial Coulombic mass to a radial Coulombic poten- 
tial. Azimuthal symmetry is assumed and spectral signatures of various 
z-dependent interaction potentials are reported. 
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1 Introduction 

The von Roos Hamiltonian [1] is known to describe quantum particles with 
position-dependent-mass (PDM), M (r) — m m(r). Over the last few decades, 
the position-dependent-mass Hamiltonians have inspired research attention [2- 
32] because of their applicability in the study of many-body problem, semicon- 
ductors, quantum dots, quantum liquids, etc. The kinetic energy operator in 
the von Roos Hamiltonian (with m Q = h = 1 units) 



H = - 1 - 

4 



m(r)~<Vm(f) p -Vm(r) a +m{r) a V m (r) f3 -V m (f)~< +V(r), (1) 



admits an ordering ambiguity manifested by the non-uniqueness representation 
of the kinetic energy operator. Which would, in effect, introduce a profile change 
in the effective potential as the values of the parameters a, /?, and 7 change [cf., 
e.g., 25-29]. Here, a, (3, and 7 are called the von Roos ordering ambiguity 
parameters satisfying the von Roos constraint a + /3 + 7 = — 1. Nevertheless, 
an interesting and comprehensive background on the issue of consistency and 
usefulness of the position-dependent mass Schrodinger equation was given by 
Levy-Leblond [32] . Therein, his work is devoted to sustaining and strengthening 
the conclusions that not only the use of position-dependent mass gives correct 
approximation, but its also a conceptually consistent approach. 

It is however constructive to mention that the continuity conditions at the 
abrupt heterojunction between two crystals implied that a = 7, otherwise for 
a^7 the wavefunctions vanish at the boundaries and the heterojunction plays 
the role of impenetrable barrier (cf., e.g., Mustafa and Mazharimousavi [10] 
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and Koc et al. in [28]). Eliminating in the process, the Gora's and Williams' 
(f3 = 7 = 0, a = —1), and Li's and Kuhn's ((3 = 7 = —1/2, a = 0) known 
parametric sets. Moreover, Dutra's and Almeida's [9] reliability test classifies 
the parametric set of Ben Daniel and Duke (a = 7 = 0, /3 = — 1) as a set to-be- 
discarded for it yields imaginary eigenvalues. This would leave us with Zhu's 
and Kroemer's (0 = 7 = —1/2, /3 = 0) and Mustafa's and Mazharimousavi's 
(a = 7 = — 1/4, /3 = —1/2) ordering ambiguity parameters that are classified 
as "good" parametric sets, so to speak. Nevertheless, we have observed (cf., e. 
g., [29]) that the physical and/or mathematical admissibility of a given ambi- 
guity parametric set depends also on the form of the position-dependent-mass 
and/or the form of the interaction potential. In the forthcoming methodical 
proposal, we shall work with the ambiguity parameters as they are without any 
classification as to which set is "good" or "to-be-discarded". 

Very recently, Mustafa [31] has considered the von Roos Hamiltonian (1) 
using cylindrical coordinates. Therein, we sought some manifestly feasible sep- 
arability through the suggestion that the position-dependent-mass is only radial- 
dependent (i.e., m(r) = m M (p,ip, z) = M(p,ip,z) = M (p) = l/p 2 ), where 
azimuthal symmetrization is granted through a proper assumption of the inter- 
action potential. The spectral signatures of different z-depcnclent interaction 
potential settings on the radial Coulombic and radial harmonic oscillator in- 
teraction potentials' spectra are reported for impenetrable walls at z = and 
z = L, for a Morse, for a non-Hcrmitian PT-symmetrized Scarf II, and for a 
non-Hcrmitian 'PT-symmetrized Samsonov interaction models. 

In this work, we offer a parallel azimuthal symmetrization along with a more 
general (though still only radial-dependent) power-law-type position-dependent- 
mass (i.e., M(p,<p,z) = M (p) = bp 2v+1 /2). Obviously, a v = -3/2 and b = 2 
yield M (p) ~ l/p 2 which is, under the current forthcoming settings, a special 
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case of M (p) = bp 2v+1 /2 that shall not be repeated here. Instead, we shall 
use v = — 1 and v = 1/2 that yield quantum particles endowed with position- 
dependent masses of a Coulombic-type, M (p) = bp~ 1 /2, and a harmonic os- 
cillator type, M (p) = bp 2 /2, respectively. To the best of our knowledge, such 
position-dependent mass settings have not been considered elsewhere. 

To make this work self-contained, we recollect (in section 2) the most rele- 
vant and vital relations (namely, equations (2)-(5) below) that have been readily 
reported in [31] for cylindrical coordinates separability and exact solvability of 
the PDM-Hamiltonian (1). In the same section, we discuss the separability in 
the framework of a manifestly radial power-law position-dependent mass and 
contemplate on the feasible separabilities. In section 3, we consider two par- 
ticular radial mass settings; a harmonic oscillator type, M (p, tp,z) = M (p) = 
g (p) = bp 2 /2, and a Coulombic type, M (p, tp, z) = M (p) = g (p) = bp~ 1 /2. We 
subject the radial harmonic oscillator type mass to move under the influence 
of a radial harmonic oscillator potential field V (p) = a 2 p 2 /4 and the radial 
Coulombic mass to a radial Coulombic potential V (p) = — 2 A/ p. The spectral 
signatures of (i) two impenetrable walls at z — and z = L provided by the 
potential well V (z) = for < z < L and V (z) — oo elsewhere, (ii) a Morse 
type [31] interaction V (z) = D (e~ 2ez — 2e~ £Z ) ; D > 0, and (iii) a trigonomet- 
ric Rosen-Morse [32] potential V (z) = U cot 2 (nz/d) :z£ [0, d\, are reported in 
the same section. Our concluding remarks are given in section 4. 

2 Cylindrical coordinates and radial power-law 
PDM framework 

Following our recent work [31] on cylindrical coordinates separability and ex- 
act solvability of the PDM-Hamiltonian (1), we again consider the position- 
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dependent-mass and the interaction potential to take the forms m(r) = M (p, p, z) — 
g (p) f (ip) k (z) and V (r) = V (p, tp, z), respectively. We have shown (see Mustafa 
[31] for more details on this issue) that the corresponding PDM-Schrodinger 
equation [H — E] ^ (p, ip, z) — with 

* (p, <p,z)=R (p) $ {ip) Z(z); p e (0, oo) , Lp e (0, 2tt) , z E (-oo, oo) , (2) 

would imply 

= 2g(p)f(<p)k(z)[E-V(p,<p,z)] 
\R"{p) (g'{p) 1\ R'jp) 
[R(p) \g(p) p) R(p) 

, C (9'(p)\ 2 W + i) ( 1(pI , g" (p) \ 

2\g{p)) 2 \pg{p) + g(p) ) 

\ Z"(Z) k^z)Z-{z)_ ( (k^V W + l)k" (z) ' 
Z(z) k(z) Z(z) 2\k(z) J 2 k(z) 

*"(<p) c(f'(v)\ 2 W + l)f"(p) } 

$M f(tp) + 2 \f{sp)) 2 f(p)\ { ) 

where 

C = a(a-l)+ 7 ( 7 -l)-/?(/3 + l). (4) 

To facilitate and ease separability, we have suggested that the interaction po- 
tential satisfies an obviously " manifested-by-equation (3)" general identity of 
the form 

2MV (p, p, z) = 2g (p) f (p) k (z) V (p, <p,z) = V (p) + V (z) + ±V (p) . (5) 

Hereby, we may remind the reader that in [31} we have used g (p) = 1/p 2 along 
with / (p) = 1 = k (z) as one of the options that secured separability of the 
problem at hand. 



+ — 

p 2 
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In the search for a more general recipe, however, we choose to eliminate 
the first-order derivatives Z' (z), ((f), and R' (p). At this point, the elimi- 
nation of the first-order derivatives of Z (z) and 3> (tp) is achieved through the 
substitutions 

Z(z)= VW)Z(z) and$(^) = V/M^M, (6) 



to imply that 



Z"(z) k'{z)Z'{z) 3{k'(z)\ 2 Ik" {z) Z"(z) 



Z{z) k{z) Z(z) A\k{z)J 2 k{z) Z(z) 

and 

$" (p) r (p) *' (p) _ 3 / /' M \ 2 , i /" m , $" m 



(7) 



$fo>) /M $M 4 V/M; + 2 + $(^) (8) 

Whereas, the elimination of the first-order derivative of R (p) may be sought 
through the substitutions 



R (p) = p u U (p) and g (p) = ^p 2v+1 ;v, b e R, (9) 



(with the restriction that b is a non-zero constant to avoid triviality) to imply 
that 

R"{p) _ (?M i\ #W = U^(P± tjOj+i) (m) 
R( P ) \g(p) p) R{ P ) U(p) p* ■ [ > 

It should be noted here that the choice of g (p) = bp 2v+1 /2 in (9) is manifestly 
mandated by the elimination of the first-order derivative of U (p). 
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Under such settings, equation (3) would read 



= bp 2v+1 f{p)k{z)E 



U"{p) (2w + l) 2 [C- P-\\-2v(v + \) 



1 



U{ P ) 

Z"{z) 
Z(z) ~" 

(g) 



2p 2 



V(p) 



(2C-3) f k> (z) 
k{z) 



P k (z) 
2 k{z) 



V(z) 



+ 



(2C-3) f f (<p) 
f(<P) 



2 /(¥>) 



(11) 



Hereby, it should mentioned that the case where t; = —3/2, 6 = 2 (i.e., 5 (p) = 
1/p 2 ) and fc (0) = 1 = / ((f) is the case we have considered in [31]. It is now just 
a special case of the current, though more general, recipe g (p) = bp 2v+1 /2; v, b e 
R. The results and examples reported therein are recoverable and hold true for 
the current work, therefore. 

Yet, an obvious manifestation of the energy term bp 2v+l f (ip) k (z) E towards 
separability of (11) is that, in addition to the three feasible separable cases 
f{<P) = 1 = k ( z ), k(z) = 1 = g(p), and / (ip) = 1 = g(p) (reported in [31]), 
one finds two more feasibly separable cases. They are, for v = —3/2, k (z) = 1, 
/ (ip) 7^ 1 (which, in turn, would break azimuthal symmetry) and for v — —3/2, 
k(z) 7^ 1, / (ip) = 1. Therefore, the separability of (11) may be facilitated 
by the forms of the position-dependent mass and of the interaction potential 
V(p,ip,z). 

To secure azimuthal symmetrization of the problem at hand we substitute 
V (ip) — and / (ip) = 1. Moreover, we choose k (z) = 1 to imply that 



k 2 ■ k 2 



-m 



m 



0,1,2, 



(12) 



d 2 z +V(z) Z(z) = k 2 z Z(z), 



(13) 
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and 



-dl + ^4^ + V{p)-b P 2v+1 E 



U (p) = -k 2 z U (p) . 



(14) 



Where to is the magnetic quantum number and 



f 



9 1 (2v + iy[(-p- ll 

' v (v + 1) + m 2 + - - ^ i. 



(15) 



is an irrational magnetic quantum number. Hereby, it obvious that the substi- 
tutions of v = —3/2 and b = 2 in (14) would inspire a re-scale of the form 



-3/2 



+ 2E = (to 2 + 3) - 2 (C - /?) + 2E, 



(16) 



so that our results in [31] for the radial Coulombic V (p) = —2/p and the radial 
harmonic oscillator V (p) = a 2 p 2 /4 (equations (23) and (24) in [31], respectively) 
are safely reproduced. Therefore, the corresponding spectral signatures of V (z) 
interaction potentials for v = —3/2 and 6 = 2 (i.e, g (p) = p~ 2 ) are reported 
therein [31] and shall not be repeated here. 



3 Two particular radial settings; g (p) = bp 2 / 2 

and g (p) = bp^/2 

In this section, we consider the position-dependent mass M (p, ip,z) = g (p) 
to indulge a radial harmonic oscillator g(p) = bp 2 /2 (i.e., v — 1/2) and the 
radial Coulombic g (p) = bp^ 1 /2 (i.e., v = —1) forms. For the sake of keeping 
this work simple and instructive, we shall consider the radial harmonic oscilla- 
tor g (p) = bp 2 / 2 accompanied by a radial harmonic oscillator type interaction 
V (p) = a 2 p 2 /A and the radial Coulombic g (p) = bp^ 1 /2 accompanied by a 
radial Coulombic V (p) = — 2 A/ p. We shall moreover report the spectral signa- 
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tures of different V (z) potentials on the overall spectrum. 

3.1 The radial harmonic oscillator g (p) = bp 2 /2 

The choice of v = 1/2 along with V (p) = a 2 p 2 /4: would imply that equation 
(14) reads 



U(p) = -klU(ph (17) 



and (15), in turn, yields 



'■1/2 



= v / (m 2 + 3)-2(C- / 3). 



(18) 



Obviously, Eq.(17) has exact eigenvalues in the form 



k 2 z = -^(a 2 -4bE) 



2n p + 



+ 1 



(19) 



and implies that 



E= a --^ 
46 Ab 



2n p + ^(m 2 + 3) - 2 (C - /3) + 1 



(20) 



We observe that an auxiliary constraint 



(C - P) - a (a - 1) + 7 ( 7 - 1) - (3 ((3 + 2) < (m 2 + 3) /2 (21) 



on the ambiguity parameters is manifested here by the requirement that E £ K. 



3.1.1 Spectral signatures of some V (z) potentials on the radial har- 
monic oscillator spectrum 

Recollect [31] that if our PDM-particle is trapped to move between two impen- 
etrable walls at z = and z = L under the influence of a 



; < z < L 
1 (.--)=<! , (22) 

oo ; elsewhere 



one would find that K z — n z Tr/L , n z = 1, 2, 3, • • • (see [31] for more details on 
this issue). This would, in effect, give the spectral signature of V (z) of (22) on 
the overall spectrum 



4b 4b 



(n z n/L) 2 



2n p + ^(m 2 + 3) - 2 (C - J) + 1 



(23) 



for a PDM particle of M (p, ip,z) = M (p) — bp 2 /2 moving in a potential of the 
form 

a 2 I ; < z < L 
V(j>,<p,z) = -+\ . (24) 

I oo ; elsewhere 

Next, let us subject this PDM particle to move in a Morse type [31] inter- 
action V(z) = D (e- 2ez - 2er cz ) ; D > 0. In this case 



^--n z - \ j , h z =0,1,2,3,--- (25) 



Therefore, a PDM quantum particle endowed with M (p, <p,z) = M (p) — bp 2 /2 
and subjected to an interaction potential of the form 

V{p,<p,z) = £ + -^ (e- 2 " - 2e"«) ;D>0 (26) 
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would admit exact energy eigenvalues given by 



E = 



1 

46 ' 46 



e-n. 



2n p + ^(m 2 + 3) -2(C-/8) + l 



(27) 



Now, let M (p, <p, 0) = M (p) = 6p 2 /2 move under the influence of a trigono- 
metric Rosen-Morse potential V (z) = U cot 2 (nz/d) ; z e [0, d] . Where ?7 and 
d are two positive parameters. In this case, 



a 2 U 

V{p,<p,z) = — + -^ cot 2 (irz/d) ;ze[0,rf], 



(28) 



*, 2 = ^[Cd + n,*] 2 -tfo;C=£ 1 + ^1 + ^ 



(29) 



(sec Ma ct al [32] for more details, notice that one should consider 2p = h = 1 
of Ma as proper parametric mapping into our settings) and 



E : 



46 46 



[Cd + n z irf jd 2 - U t 



-1 2 



2n p + V(™ 2 + 3) - 2 (C - /3) + 1 



(30) 



3.2 The Radial Coulombic g (p) = bp~ 1 /2 



Now consider the PDM-particle to have a radial Coulombic-type mass of the 
form M(p,<p,z) = M (p) = 6p~ 1 /2, (i.e., u = —1) and subjected to move in a 
radial Coulombic potential V (p) = — 2 A/ p. In this case, 



f (p, ¥>,*) = -y + £^(*)> 



(31) 



and equation (14) yields 



"3 + 



.i-1/4 2(6£; + ,4 



£/(p) = -fc 2 *7(p); 6 = 6/2, (32) 
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with 



= v/(™ 2 + 3/4)-(C-/?) A 



(33) 



and 



Which, in turn, results 

E = ±j (n p + + 3/4) - (C - 0) /2 + l) - |, (35) 

with the auxiliary constraint 

(C-/3)=a(a-l)+7(7-l)-/3(/8 + 2)< (2m 2 + 3/2) , (36) 

on the ambiguity parameters that secures the reality of E. Nevertheless, two 
branches of energies are obviously obtained. Moreover, the spectral signature 
of k z on the overall spectrum is obtained through the solution of equation (13). 

3.2.1 Spectral signatures of some V (z) potentials on the radial Coulom- 
bic spectrum 

If we subject our radial Coulombic PDM-particle M (p, tp,z) = M (p) = bp^ 1 /2 
to move in V (p, f, z) = —A/b + pV (z) /b, where V (z) is given by (22), it will 
admit exact energies of the form 

E = ±5f (n p + v/(m 2 + 3/4) - (C - J) /2 + l) - i; n z = 1, 2, 3, • • • . (37) 
bL v ' b 

Moreover, if this PDM-particle is subjected to move in a Morse type [31] 



k z = ±- 



bE + A 
+ I-i 



+ 1 



(34) 
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interaction V (z) = D (e~ 2 " - 2e~ £Z ) ; D > 0. In this case, 



V (p, <p, z) = -j + \D ( e - 2ez - 2e-«) ; .D > 0, 



(38) 



and the exact energies are of the form 



E = ±^ 



- n, - ij (», + V(™ 2 + 3/4)-(C-ffl/2 + l) - 4, (39) 



Whcrcri2 = 0, 1,2,3, Obviously, the condition ^v^D/e — n z — > is 
manifested here and ought to be enforced, otherwise complex pairs of energy 
eigenvalues are obtained in the process. 

Next, let M(p,ip,z) — M (p) = bp~ 1 /2 move under the influence of a 
trigonometric Rosen- Morse potential V (z) = U Q cot 2 (nz/d) ;z£ [0,d]. Then, 



A p 

V(p, <p,z) = -— + jU cot 2 (nz/d) ;z <= [0,d] , 



(40) 



and 



E 



J[Cd + h z TT} 2 -U d 2 ( 

= ±- f 

bd v 



+ v/(™ 2 + 3/4)-(C-/3) /2 + l) 



T- (41) 





4 Concluding remarks 

We have recollected the most relevant and vital relations (equations (2)-(5) 
above) that have been readily reported by Mustafa [31] for cylindrical coor- 
dinates separability of the PDM-Hamiltonian in (1), where the PDM-setting 
was considered in the form M (p, <p, z) = g (p) f (ip) k (z) = g (p) = 1/p 2 under 
azimuthally symmetric settings. 

In this work, however, we offered a more general power-law radial position- 
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dependent mass recipe M(p,ip,z) = g (p) = bp 2v+1 /2; v, b e M, within which 
M (p, (p,z) = g (p) = 1/p 2 of [31] represents a special case (the results and ex- 
amples reported therein hold true and yet document additional examples on 
the applicability of the current methodical proposal, therefore). Moreover, the 
structure of the position-dependent energy term bp 2v+1 f (<p) k (z) E in (11) sug- 
gests that there are five feasible cases towards separability; (i) / (ip) = 1 = k (z), 
(ii) k(z) = l=g (p), (iii) f( !P ) = \ = g (p), (iv) v = -3/2, k (z) ± 1, / (cp) = 1, 
and (v) v = —3/2, k (z) — 1, / (tp) 1 (which would break azimuthal sym- 
metry, of course). Therefore, the separability of (3) may be facilitated by the 
forms of the position-dependent mass and the interaction potential V (p, ip, z). 
These are not the only cases to secure separability of (3), so to speak. 

We have considered two particular mass settings; a radial harmonic oscil- 
lator type, M (p, <p,z) — M (p) = g (p) = bp 2 /2, and a radial Coulombic type, 
M (p, <p,z) — M (p) — g(p) = bp~ 1 /2. We have observed that for the Coulombic 
case two branches of energies are obtained, each of which is a "mirror-reflection" 
of the other about the zero-energy axis. Moreover, when we subjected the radial 
harmonic oscillator mass to a radial harmonic oscillator potential V (p) = a 2 p 2 /4 
and the radial Coulombic mass to a radial Coulombic potential V (p) = —2 A/p, 
only constant shifts in the energies where observed (i.e., a shift (a 2 /4&) for the 
radial harmonic oscillator mass and ^— A/bj for the radial Coulombic mass, 
documented in (20) and (35), respectively). That is, the radial interaction po- 
tentials V (p) considered for the two over simplified examples here provided no 
quantization recipe at all (i.e., they have only introduced constant shifts to the 
energies but not discrete quantum energy shifts). This is because the form of 
the general interaction potential V (p, ip, z) we have adopted in (5). 

Yet, auxiliary constraints on the ambiguity parameters (see (21) for the 
harmonic oscillator and (36) for the Coulombic) are observed mandatory to 
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secure the reality of E. Hereby, if m = is considered in (21) and (36) as 
a reference test, then one would observe that only the Gora's and Williams' 
ambiguity parametric set (ft = 7 = 0, a = — 1) fails to provide real energies 
(i.e., \/3/2 — (C — /?) € C) . We contemplate that more auxiliary constraints 
on the ambiguity parameters should be anticipated for different, though ex- 
actly solvable, power-law type radial masses (within our methodical proposal, 
of course). Furthermore, the spectral signatures of different V (z) interactions 
on the overall spectrum are also reported. Namely, the spectral signatures of 
(i) two impenetrable walls at z = and z — L provided by the potential well 
V (z) — for < z < L and V (z) = 00 elsewhere, (ii) a Morse type [31] interac- 
tion V (z) — D (e~ 2ez — 2e~ ez ) ; D > 0, and (hi) a trigonometric Rosen-Morse 
[32] potential V (z) = U D cot 2 (irz/d) ; z e [0, d}. 
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